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We present measurements of the branching fraction and fraction of longitudinal polarization for the 
decay B+ -> K*°K*+ with a sample of 467 ± 5 million BB pairs collected with the BABAR detector 
at the PEP-II asymmetric-energy e^e~ collider at the SLAC National Accelerator Laboratory. We 
obtain the branching fraction B{B+ K''°K*+) = (1.2 ± 0.5 ± 0.1) x lO"*^ with a significance of 
3.7 standard deviations including systematic uncertainties. We measure the fraction of longitudinal 
polarization /l — 0.75l^Q'2g ± 0.03. The first error quoted is statistical and the second is systematic. 

PACS numbers: 13.25.Hw, ll.30.Er, 12.15.Hh 



The study of the branching fractions and angu- 
lar distributions of B meson decays to hadronic final 
states without a charm quark probes the dynamics of 
both the weak and strong interactions, and plays an 
important role in understanding CP violation in the 
quark sector. Improved experimental measurements of 



these charmless decays, combined with theoretical de- 
velopments, can provide significant constraints on the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix parame- 
ters [1\ and uncover evidence for physics beyond the Stan- 
dard Model 

QCD factorization models predict the fraction of Ion- 
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gitudinal polarization Jl of the decay of the B me- 
son to two vector particles {VV) to be ~ 0.9 for both 
tree- and loop-dominated (penguin) decays However, 
measurements of the penguin VV decays (f)K*^ 
and — *■ (jiK*^ give /l approximately 0.5 while 
Jl~ 0.8lto'i2 ^ 0-06 has been measured for the decay 
S° x*^K*° 0. Several attempts to understand the 
values of /l within or beyond the Standard Model have 
been made 0. Further information about decays re- 
lated by SU{3) symmetry may provide insights into this 
polarization discrepancy and test possible modifications 
to factorization models, such as penguin annihilation or 
rescattering [§]. 

The decay B+ k*^K*+ occurs through both elec- 
troweak and gluonic b ^ d penguin loops, as shown in 
Fig. [TJ Its branching fraction is expected to be of the 
same order as B° — > K*°K*°, with Beneke, Rohrer and 
Yang predicting (O.^toltoi) x IQ-^, while Cheng and 
Yang ;.3:] quote (0.6 ± 0.1 ± 0.3) x 10_-^ both based on 
QCD factorization. The B° K*°K*° branching frac- 
tion has been measured to be (1.28;gjg±0.11) x 10"^ Q, 
where the first error is statistical and the second system- 
atic, while an upper limit at the 90% confidence level 
(C.L.) of 2.0 x 10^^ has been recently placed on the 
K*~K*+ branching fraction The current ex- 
perimental upper limit on the i?+ ~^ K*'^K*~^ branching 
fraction at the 90% C.L. is 71(48) x lO'^ 0, assuming 
a fully longitudinally (transversely) polarized system. 




FIG. 1: The electroweak (left) and gluonic (right) b ^ d 
penguin loop diagrams for K'^^K*'^ . 



We report on a search for the decay mode 
K*°K*+, where K* refers to the if* (892) resonance, 
with consideration of nonresonant backgrounds [Il[ . The 
analysis is based on a data sample of 467 ± 5 million 
BB pairs, corresponding to an integrated luminosity 
of 426 fb~^, collected with the BaBAR detector at the 
PEP-II asymmetric-energy e+e^ collider operated at 
the SLAC National Accelerator Laboratory. The e"*"e~ 
center-of-mass (cm.) energy is ^/s = 10.58 GeV, cor- 
responding to the T{4:S) resonance mass (on-resonance 
data). In addition, 44.4 fb"^ of data collected 40 MeV be- 
low the T{4:S) resonance (off-resonance data) are used for 
background studies. We assume equal production rates 
of B+B- and mesons. 

The BaBAR detector is described in detail in Ref. [l3] . 
Charged particles are reconstructed as tracks with a five- 
layer silicon vertex detector and a 40-layer drift cham- 
ber inside a 1.5 T solenoidal magnet. An electromagnetic 



calorimeter (EMC) comprised of 6580 CsI(Tl) crystals is 
used to identify electrons and photons. A ring-imaging 
Cherenkov detector (DIRC) is used to identify charged 
hadrons and to provide additional electron identification 
information. The average K-tt separation in the DIRC 
varies from 12 cr at a laboratory momentum of 1.5 GeV/c 
to 2.5 a at 4.5 GeV/c. Muons are identified by an instru- 
mented magnetic-flux return (IFR). 

The B+ j^*Oj^*+ candidates are reconstructed 
through the decays of X*" K^tt+ and K*+ A'°7r+ 
or K*+ K+TT^, with if" ^ tt+tt- and 7r° jj. The 
differential decay rate, after integrating over the angle 
between the decay planes of the vector mesons, for which 
the acceptance is uniform, is 
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T d cos 6*1 d cos 6*2 

— sin^ 9i sin^ 62 + cos^ Oi cos^ 62, 



(1) 



where 9i and 6*2 are the helicity angles of the if *+ and 
if *°, defined as the angle between the daughter kaon (if" 
or if^) momentum and the direction opposite to the B 
meson in the if* rest frame [isj . 

The charged particles from the if* decays are required 
to have a transverse momentum relative to the beam axis 
greater than 0.05GeV/c. The particles are identified as 
either charged pions or kaons by measurement of the en- 
ergy loss in the tracking devices, the number of photons 
recorded by the DIRC and the corresponding Cherenkov 
angle. These measurements arc combined with additional 
information from the EMC and IFR detectors, where ap- 
propriate, to reject electrons, muons, and protons. 

The if g candidates are required to have a mass within 
0.01 GeV/c^ of the nominal if" mass 14], a decay vertex 
separated from the B meson decay vertex by at least 
twenty times the uncertainty in the measurement of the 
separation of the vertex positions, a flight distance in the 
direction transverse to the beam axis of at least 0.3 cm, 
and the cosine of the angle between the line joining the 
B and if" decay vertices and the if" momentum greater 
than 0.999. 

In the laboratory frame, the energy of each photon 
from the tt" candidate must be greater than 0.04 GeV, 
the energy of the tt" must exceed 0.25 GeV, and the re- 
constructed tt" invariant mass is required to be in the 
range 0.12 < m^^ < 0.15 GeV/ c^. After selection, the tt" 
candidate's mass is constrained to its nominal value 14| . 

We require the invariant mass of the if* candidates 
to satisfy 0.792 < m^^ < 0.992 GeV/c^. A B meson 
candidate is formed from the if*" and A'*+ candidates, 
with the condition that the if* candidates originate from 
the interaction region. 

The B meson candidates are characterized kincmat- 
ically by the energy difference AE — Eg — \/s/2 
and the beam energy-substituted mass = 
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[(s/2 + p, •pB)VKf ~p|]'/', where {E„pi) and 
{Eb,Pb) are the four-momenta of the T(4S') and B 
meson candidate in the laboratory frame, respectively, 
and the asterisk denotes the cm. frame. The total event 
sample is taken from the region —0.10 < AE < 0.15 GeV 
and 5.25 < toes < 5.29 GeV/ c^. The asymmetric 
AE criterion is applied to remove backgrounds from 
B to charm decays that occur in the negative AE 
region. Events outside the region \AE\ < 0.07 GeV and 
5.27 < Toes < 5.29 GeV/ are used to characterize the 
background. 

We suppress the background from decays to charmed 
states by forming the invariant mass uid from combina- 
tions of three out of the four daughter particles' four- 
momenta. The event is rejected if 1.845 < mu < 
1.895 GeV/ and the charge and particle type of the 
tracks are consistent with a known decay from a D me- 
son [13] ■ Backgrounds from B (jjK* are reduced by 
assigning the kaon mass to the pion candidate and re- 
jecting the event if the invariant mass of the two charged 
tracks is between 1.00 and 1.04GeV/c^. Finally, to re- 
duce the continuum background and to avoid the region 
where the reconstruction efficiency falls off rapidly for low 
momentum tracks, we require the cosine of the helicity 
angle of the K* candidates to satisfy cos 9 < 0.98. 

To reject the dominant background consisting of light- 
quark qq (q = u,d,s,c) continuum events, we require 
I cos 6't I < 0.8, where 9t is the angle, in the cm. frame, 
between the thrust axis [l3] of the B meson and that 
formed from the other tracks and neutral clusters in the 
event. Signal events have a flat distribution in |cos0t|j 
while continuum events peak at 1. 

We use Monte Carlo (MC) simulations of the signal 
decay to estimate the number of signal candidates per 
event. After the application of the selection criteria, the 
average number of signal candidates per event is 1.06 
(1.02) for fully longitudinally (transversely) polarized de- 
cays with no tt" in the final state and 1.15 (1-07) for de- 
cays with one 7r° in the final state. The candidate with 
the smallest fitted decay vertex is chosen. MC sim- 
ulations show that up to 5.1% (1.7%) of longitudinally 
(transversely) polarized signal events with no 7r° are mis- 
reconstructed, with one or more tracks originating from 
the other B meson in the event. In the case of signal 
events with one 7r°, the fraction of misreconstructed can- 
didates is 8.8% (2.8%) for longitudinally (transversely) 
polarized signal events. 

A neural net discriminant J\f is used in the maximum- 
likelihood (ML) fit, constructed from a combination of six 
variables calculated in the cm. frame: the polar angles of 
the B meson momentum vector and the B meson thrust 
axis with respect to the beam axis, the angle between 
the B meson thrust axis and the thrust axis of the rest 
of the event, the ratio of the second- and zeroth-order 
momentum-weighted Legendre polynomial moments of 
the energy flow around the B meson thrust axis |16l] , the 



flavor of the other B meson as reported by a multivariate 
tagging algorithm [17|, and the boost-corrected proper- 



time difference between the decays of the two B mesons 
divided by its variance. The second B meson is formed by 
creating a vertex from the remaining tracks that are con- 
sistent with originating from the interaction region. The 
discriminant is trained using MC for signal, and qq con- 
tinuum MC, off-resonance data and on-resonance data 
outside the signal region for the background. 

An extended unbinncd ML fit is used to extract the sig- 
nal yield and polarization simultaneously for each mode. 
The extended likelihood function is 



N 



c 



Nl 



exp 



E". n 



(2) 



We define the likelihood £i for each event candidate i as 
the sum of njVj{xi; dj) over three hypotheses j: sig- 
nal (including misreconstructed signal candidates), qq 
background and BB backgrounds as discussed below. 
Vj{xi; dj) is the product of the probability density func- 
tions (PDFs) for hypothesis j evaluated for the i-th 
event's measured variables the yield for hypoth- 

esis j , and N is the total number of events in the sam- 
ple. The quantities dj represent parameters to describe 
the expected distributions of the measured variables for 
each hypothesis j . Each discriminating variable Xi in the 
likelihood function is modeled with a PDF, where the 
parameters dj are extracted from MC simulation, off- 
resonance data, or (toeSj AE) sideband data. The seven 
variables Xi used in the fit are toes, AE, Af, and the in- 
variant masses and cosines of the helicity angle of the two 
K* candidates. Since the linear correlations among the 
fitted input variables are found to be on average about 
1%, with a maximum of 5%, we take each Vj to be the 
product of the PDFs for the separate variables. 

For the signal, we use relativistic Breit-Wigner func- 
tions for the K*^ and K*^ invariant masses and a sum of 
two Gaussians for toes and AE. The longitudinal (trans- 
verse) helicity angle distributions are described with a 
cos^ 9 (sin^ 9) function corrected for changes in efficiency 
as a function of helicity angle. The correction also ac- 
counts for the reduction in efficiency at a helicity near 
0.78 introduced indirectly by the criteria used to veto 
D mesons. The continuum toes shape is described by 
the function xVl — x"^ exp[— £(1 — cc^)] with x = m-Es/Eg 
and £ a free parameter |18j , while second- and third-order 
polynomials are used for AE and the helicity angles, re- 
spectively. The continuum invariant mass distributions 
contain peaks due to real K* candidates; we model the 
peaking mass component using the parameters extracted 
from the fit to the MC signal invariant mass distributions 
together with a second-order polynomial to represent the 
nonpeaking component. The BB backgrounds use the 
same toes function as the continuum and an empirical 
nonparametric function [l9l| for all other variables. The 
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neural net distributions are modeled using the empirical 
nonparametric function for all hypotheses. 

BB backgrounds that remain after the event selection 
criteria have been applied are identified and modeled us- 
ing MC simulation [20|. There are no significant charm- 
less BB backgrounds. However, decays from higher mass 
fC*°(1430) states are not fully simulated due to their 
uncertain cross-section and resonance structure and we 
treat these as an explicit systematic uncertainty later. 
The charm BB backgrounds are effectively suppressed 
by applying the veto on the D meson mass described 
above. The remaining charm BB background events are 
mostly single candidates formed from the decay products 
of a D, D* or D*^, together with another track from the 
event. The polarization and branching fractions of these 
backgrounds are uncertain and so we fix the BB back- 
ground yield in the fit and then vary the yield by ±100% 
as a systematic cross-check. 

We fit for the branching fraction B and simultane- 
ously and exploit the fact that B is less correlated with Jl 
than is either the yield or efficiency taken separately. The 
continuum background PDF parameters that are allowed 
to vary are ^ for tties, the slope of AE, and the polyno- 
mial coefficients and normalizations describing the mass 
and helicity angle distributions. We validate the fitting 
procedure and obtain the sizes of potential biases on the 
fit results by applying the fit to ensembles of simulated 
experiments using the extracted fitted yields from data. 
The qq component is drawn from the PDF, and the signal 
and BB background events are randomly sampled from 
the fully simulated MC samples. Any observed fit bias is 
subtracted from the fitted yield. 

The total event sample consists of 1381 and 3201 
events for B+ k*Ox*+ with K*+ -> K'^t:+ and 
K*+ K+TT°, respectively. The results of the ML fits 
are summarized in Table ID We compute the branch- 
ing fractions B by dividing the bias-corrected yield by 
the number of BB pairs, the reconstruction efficiency e 
given the fitted /l, and the secondary branching frac- 
tions, which we take to be 2/3 for B{K*" K-'k+) 
and B{K*+ K°Tr+), 1/3 for B{K*+ K+Tr°), and 
0.5 X (69.20 ± 0.05)% for B{K° -> ifj? -> ir+n -). The 
significance S of the signal is defined as 5 = ■\/2Aln£, 
where A In £ is the change in log- likelihood from the max- 
imum value when the number of signal events is set to 
zero, corrected for the systematic errors by convolving 
the likelihood function with a Gaussian distribution with 
a variance equal to the total systematic error defined be- 
low. We confirm that 2Aln£ is a reliable estimate of 
the significance S by fitting ensembles of simulated ex- 
periments with background events only, using the fitted 
parameters and background yields from the data, and ob- 
serving how often the number of fitted signal events ex- 
ceeds the fitted signal yield in the data. The significance 
of the combined B"*" branching fractions is 

3.7cr, including statistical and systematic uncertainties. 



The 90% C.L. branching fraction upper limit (iJuL) is de- 
termined by combining the likelihoods from the two fits 
and integrating the total likelihood distribution (taking 
into account correlated and uncorrelated systematic un- 
certainties) as a function of the branching fraction from 

to KuL, so that /(f""^ CdB = 0.9 CdB. 

Figures [2] and [3] show the projections of the two fits 
onto TOES, Ai?, and the K*^ and K*^ masses and 
cosines of the helicity angle for the final state with zero 
and one tt*^, respectively. The candidates in the figures 
are subject to a requirement on the probability ratio 
Psig/('Psig + ^bkg), optimized to enhance the visibility 
of potential signal, where Vsig, and Pbkg are the signal 
and the total background probabilities, respectively, com- 
puted without using the variable plotted. The dip in he- 
licity at ~ 0.78 is created by the criteria used to veto the 
charm background. 

TABLE L Summary of results for the fitted yields, fit bi- 
ases, average reconstruction efficiencies e for the fitted /l, 
sub-branching fractions W Bi , longitudinal polarization /_l , 
branching fraction B {B^ — > K*'^K*^), B significance 5, and 
90% C.L. upper limit Bul- The first error is statistical and 
the second, if given, is systematic. 



Final State A" 7r+ Kl 7r+ A" 7r+ A+ tt' 



Yields (events): 






Total 


1381 


3201 


Signal 


6.9t3-5 


13.9t^:| 


qq bkg. 


1365 ± 37 


3169 ± 57 


BB bkg. (fixed) 


10 


19 


ML Fit Biases 


-0.12 


0.08 


Efficiencies and B: 






e(%) 


11.44 ±0.08 


7.40 ± 0.08 


ns,(%) 


15.37 


22.22 




0.72lJ!;3^ ± 0.03 0.79lo;36 ± 0.03 


B (xlO-^) 


0.85ini±0.11 1.80lJ;85 ±0.16 


B Significance S (a) 


2.28 


2.18 


Combined Results: 








U. 1 O_o.26 


±0.03 


B (xlO"'') 


1.2 ±0.5 ±0.1 



B Significance S (cr) 3.7 
guL (xlO-*^) 2^ 



The systematic uncertainties on the branching frac- 
tions are summarized in Table [TTl The major uncer- 
tainty is the unknown background from the decays B'^ — > 
if*+i?*0(1430) and B+ K*+ {lASQ) . We take 

the central value of the branching fraction of B+ 
i?*0(1430)A:+ [in as an estimate of the A:*(1430) K* 
branching fraction. We use the LASS parameteriza- 
tion for the AT* (1430) lineshape, which consists of the 
AT* (1430) resonance together with an effective-range non- 
resonant component |22l |. and assume that interference 
effects between the K* and the spin-0 final states (non- 
resonant and AT* (1430)) integrate to zero as the accep- 
tance of the detector and analysis is almost uniform. This 
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FIG. 2: Projections for B+ ~* K*°{-^ K--k+)K*+{-^ J^'sTt+I 
of the multidimensional fit onto (a) tties; (b) Ai?; (c) K* 
mass; (d) cosine of K*^ helicity angle; (e) A'*" mass; and 
(f) cosine of K*^ helicity angle for events selected with a 
requirement on the signal-to-total likelihood probability ra- 
tio, optimized for each variable, with the plotted variable ex- 
cluded. The points with error bars show the data; the solid 
line shows signal-plus-background; the dashed line is the con- 
tinuum background; and the hatched region is the signal. 



lineshape is used to calculate the number of background 
events in the K* mass range. We estimate 0.81 and 0.77 
events in the modes without and with a tt" , respectively. 

The other errors on the branching fractions arise from 
the PDFs, fit biases and BB background yields, and effi- 
ciencies. The PDF uncertainties are calculated by vary- 
ing the PDF parameters that are held fixed in the original 
fit by their errors, taking into account correlations. The 
uncertainty from the fit bias includes its statistical un- 
certainty from the simulated experiments and half of the 
correction itself, added in quadrature. The uncertain- 
ties in PDF modeling and fit bias are additive in nature 
and affect the significance of the branching fraction re- 
sults. Multiplicative uncertainties include reconstruction 
efficiency uncertainties from tracking and particle iden- 
tification (PID), track multiplicity, MC signal efficiency 
statistics, and the number of BB pairs. The majority of 
the systematic uncertainties on cancel and the error 
is dominated by the uncertainty on the PDF parameters. 
This is calculated to be ±0.03 for both modes. 

In summary, we have seen a significant excess of events 
and have measured the branching fraction B{B'^ — > 
K*°K*+) = [1.2±0.5(stat)±0.1(syst)] x 10"^ and the 
longitudinal polarization Jl = O.l^to^Q ± 0.03. The 90% 
C.L. upper limit on the branching fraction is found to be 



FIG. 3: Projections for B+ K*° 
of the multidimensional fit onto 



0.5 

K'° helicity 

a) mEs; (b) AE; (c) K'^ 



mass; (d) cosine of K*^ helicity angle; (e) K*'^ mass; and (f) 



cosine of Tf*" helicity angle. The same projection criteria and 
legend are used as in Fig. [2l 



TABLE II: Estimated systematic errors on the branching frac- 
tion in the final fit. Error sources which are correlated and 
uncorrelated when combined from the two decays are denoted 
by C and U, respectively. 



Final State K~ 


7V+ Kl n+ 


K- TT+ K+ 7r° 


Additive errors (events): 






Fit Bias [U] 


0.08 


0.09 


Fit Parameters [U] 


0.10 


0.39 


LASS backgrounds [U] 


0.81 


0.77 


BB backgrounds [U] 


0.10 


0.60 


Total Additive (events) 


0.83 


1.06 


Multiplicative errors (%): 






Track Muhiplicity [C] 


1.0 


1.0 


MC Statistics [U] 


0.2 


0.2 


Number of BB pairs [C] 


1.1 


1.1 


PID [C] 


3.3 


3.3 


Neutrals Corrections [C] 




3.0 


Ks Corrections [C] 


1.4 




Tracking Corrections [C] 


2.4 


2.4 


Total Muhiplicative (%) 


4.2 


4.7 


Total B error (xlO"") 


0.11 


0.16 



B{B+ K*°K*+) < 2.0 X 10-^. These measurements 
are compatible with theoretical predictions. 
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